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Abstract 

The Two-Body Dirac equations of constraint theory are of special interest not only 
in view of applications for phenomenological calculations of mesonic spectra but also 
because they avoid no-go theorems about relativistic interactions. Furthermore, they 
provide a quantum mechanical description in a manifestly Lorentz invariant way using 
the concept of a multi-time wave function. In this paper, we place them into the con¬ 
text of the multi-time formalism of Dirac, Tomonaga and Schwinger for the first time. 
A general physical and mathematical framework is outlined and the mechanism which 
permits relativistic interaction is identified. The main requirement derived from the 
general framework is the existence of conserved tensor currents with a positive com¬ 
ponent which can play the role of a probability density. We analyze this question for 
a general class of Two-Body Dirac equations thoroughly and comprehensively. While 
the free Dirac current is not conserved, it is possible to find replacements. Improving 
on previous research, we achieve definite conclusions whether restrictions of the func¬ 
tion space or of the interaction terms can guarantee the positive definiteness of the 
currents - and whether such restrictions are physically adequate. The consequences 
of the results are drawn, with respect to both applied and foundational perspectives. 

Keywords: relativistic wave equations, current conservation, multi-time wave func¬ 
tions, Lorentz invariance, relativistic interactions, multi-particle Dirac equation 


1 Introduction 

In view of the many open questions in quantum field theory (QFT) of both conceptional and 
mathematical nature, most prominently the measurement problem and the ultraviolet (UV) 
divergencies [T] , it seems a worthwhile question to ask whether there exists an intermediary 
level between non-relativistic quantum mechanics (QM) and QFT: a relativistic iV-particle 
theory. 

In this regard, the development in constraint theory [2] has led to the formulation of the 
Two-Body Dirac (2BD) equations [SHS]. These are simultaneous relativistic equations for 
a 16-component wave function q(a;i, X 2 ) which relates the space-time coordinates Xi of two 
particles i = 1,2. They take the following form: 

Di'ip{xi,X2) = 0 , 

D2'4){xi,X2) = 0 , ( 1 ) 

where = Di^Q-\-Vi{x,P). Here, is the operator in the free Dirac equation = 0 

for the i-th particle in manifestly covariant form and Vi{x, P) are functions oi x = Xi — X 2 
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and the total momentum operator P. 

The 2BD equations are interesting with respect to several aspects: 

1. They are closed equations for bound state problems. 

2. They provide evolution equations for a multi-time wave function. 

3. They include a new mechanism of relativistic interaction. 

Concerning 1., equations of this kind are useful for applications in nuclear and particle 
physics, such as calculations of meson spectra (see e.g. IMS]). 

With regard to point 2, it seems clear that a manifestly covariant description of quantum 
phenomena in the Schrodinger picture requires the wave function to depend on all of 
the space-time coordinates Xi = of each particle i. Because of the many time 

coordinates, these relativistic wave functions are called multi-time wave functions. The idea 
of a multi-time wave function was proposed by Dirac at the very beginning of relativistic 
quantum theory as one of its key elements uadi] and it was subsequently developed 
further by Bloch |15) . Tomonaga |16| and Schwinger HZ]. More recently, it has become 
clear that multi-time wave functions are also of particular interest for the foundations 
of QM as they allow for an extension of realistic quantum theories to the relativistic 
domain. Examples include relativistic GRW theories |181I19| and relativistic Bohmian 
mechanics |2U[|21| . showing possible ways how the measurement problem can be avoided 
in a relativistic context. 

Concerning 3., one may wonder how the 2BD equations escape no-go theorems about 
direct relativistic interactions such as the “no interaction theorem” |22[I23| or the recent 
result |24| . The key feature is the inclusion of “potential” terms depending on arbitrary 
powers of the total momentum operator. Possibilities of this novel - and perhaps radical 
kind are left open by the above-mentioned no-go theorems which concern Hamiltonian 
evolution equations. 

The inclusion of potential terms of this kind raises several important conceptual questions 
which have not been clarified fully in the literature on the 2BD equations: 

1. How should one understand the 2BD equations mathematically if one includes many 
times and in addition leaves the quantum mechanical Hamiltonian framework? 

2. Does the general framework for a quantum mechanical theory have to be changed 
because of the fact that the potential terms contain differential operators up to infinite 
order? 

The first question poses itself because of (a) the inclusion of many time variables in a 
multi-time wave function and (b) the fact that the potential terms in the 2BD equations 
contain time derivatives of any order. However, not even —i-§i is an operator on the com¬ 
monly used Hilbert space as the elements of the latter are functions /(x) of space 

variables only. 

The second point concerns the question if the way how a quantum mechanical theory makes 
predictions is damaged in any way by the unconventional structure of the 2BD equations. 
In essence, one requires that the wave equations imply the existence of a current and a 
corresponding density which is positive definite so that it can play the role of a probabil¬ 
ity density. In non-relativistic QM, the famous Born rule states that is the correct 
probability density for the statistical distribution of particles in experiments. In the multi¬ 
time picture one seeks a conserved tensor current {xi,X 2 ) (here for two particles), i.e. 

= 0 = The question then is whether has a positive component, such 

as in free two-time Dirac theory. This is not evident anymore for the 2BD 
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equations and was e.g. the subject of |25) . 

In this paper, we approach the above-mentioned questions as follows: We first introduce 
the perhaps unfamiliar idea of a multi-time wave function and show that it does in fact arise 
naturally from Lorentz transformations of Wparticle wave functions from non-relativistic 
QM (sec. 12.ip . This includes a brief survey of how one can understand time evolution 
in this framework. We also outline the connection to a Hilbert space setting (sec 12.21) . 
Furthermore, we comment on how the 2BD equations evade the no-go theorems for Hamil¬ 
tonian theories (sec. 1^ . Then we introduce the mathematical idea underlying the 2BD 
equations which allows them to satisfy a certain necessary compatibility condition (sec. 
13.21) . We also discuss its physical consequences, in particular the necessity of the appear¬ 
ance of a certain covariantization of the spatial distance in the center of momentum frame 
a nd the compatibility of this fact with the theory of relativity. Subsequently, we introduce 
and focus on a concrete and important class of the 2BD equations for particle-antiparticle 
pairs which was first suggested by Sazdjian [1] (sec. 13.31) . This class is related to other 
forms of the 2BD equations like the one used by Crater and Van Alstine (see mm) such 
that our analysis is sufficiently general. We point out open mathematical issues (sec. Eai 
and propose a preliminary understanding (sec. 1331) . 

After this, we come to the main part of this work: the question whether an appropriate 
tensor current exists. To approach this question, we show that the free Dirac current is not 
conserved (sec. KU- Possible replacements can, however, be found. Nevertheless, there 
does not seem to exist a general argument that any of them is positive definite. We there¬ 
fore continue the analysis by posing the question whether further conditions can render the 
currents or equivalently, the associated scalar product, positive definite. A related analysis 
of Sazdjian |25| is discussed and found incomplete. The main open points are the questions 
of 1. how one should regard the further assumptions, (a) as restrictions on the admitted 
space of functions (sec. 14.21) or (b) as restrictions on the admissible potentials 14.31 as well 
as 2. whether the two options (a), (b) or a combination thereof is physically reasonable. 
These questions are answered and the consequences are drawn. Furthermore, we briefly 
point out that the class of gauge transformations for the 2BD equations is restricted as 
compared with the free case (sec. ED. The paper ends with a discussion of the results and 
their implications. 

2 Multi-time wave functions 

In this section, we provide some general background for multi-time wave functions and 
outline their connection with physics. This is necessary for a proper understanding of the 
2BD equations as evolution equations for a wave function 'ip{xi,X 2 ) with Xi = (fj,Xj). 

2.1 Derivation, domain and evolution equations 

First, we show how the idea of a multi-time wave function arises naturally from general con¬ 
siderations. Consider an V-particle system with single-time wave function ¥?(xi, ...,xisf,t) 
as a solution of a single-time wave equation such as Schrodinger’s or Breit’s equation |26| . 
The argument (xi, ..., xjv, t) of the wave function describes a configuration of N space 
points at a certain time t. If we want to Lorentz-transform the wave function, we expect 
from basic representation theory that the argument of the wave function should transform 
according to the inverse of the Lorenz transformation A. However, in order to have a 
Lorentz transformation act on the configuration, we need to regard it as physically syn¬ 
onymous to the following collection of four-vectors at equal time: ((f,xi), ..., (t,X]\f)). Now 
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we can act with A“^ as follows: (A“^(t, xi),(t,x^r)) = x'^),(t^,x^)) where 

in general all are different. In this way, we see that the Lorentz transformation of a 
single time wave function leads to the necessity to consider multi-time wave functions 

C X • • • X ^ S, {xi,...,xn) '—^ ifixi, ..■,X]\[)- (2) 

^ V 

N times 

The connection to single-time wave functions is given by 

(/?(xi, ...,XAr,t) = 'lp{t,Xi, ...,t,XN). (3) 

Because the Lorentz transformation of a simultaneous configuration yields a space-like 
configuration, the natural domain of ijj is the set of space-like configuration^, i.e. 

^ = {{xi,...,xn) ■ {xi-Xj)f,{xi-XjY <f)'ii^ j], ( 4 ) 

Our sign convention for the flat-space metric is g = diag(l, —1, —1, —1). 

Dirac’s original suggestion for multi-time evolution equations was to prescribe N simulta¬ 
neous systems of first order partial differential equations |13| : 

i-^if = Hi'ip 
utl 


= HniP, (5) 

OtN 

where the H^, k = 1,..., A^ are differential operators on an appropriate function space and 
do not contain time derivatives. Dirac’s idea was to prescribe eq. © on but it became 
clear shortly after that this is problematic for two reasons |151I16| : Firstly, the dynamics 
may exist on =5^ but not on R^-^ and secondly, the statistical meaning of the wave function 
can in general only be valid on space-like hypersurfaces. 

For a single-time wave function, one usually prescribes initial data (/?(xi, ...,X7v,0) = 
(/9o(xi,...,x^r). Using the connection between single- and multi-time wave functions ([3]) 
we translate this into: 

V’(0,Xi, ...,0,XAr) = VJo(xi, ...,XAr). ( 6 ) 

Alternatively, one can choose to prescribe initial data on a general set where S is a 
space-like hypersurface. One re-obtains a single-time wave equation 



N 


k=l 


(7) 


from eqs. ([5]) by using the chain rule together with eq. ©• 

Because of the many simultaneous equations ([5]), questions of consistency arise. For the 
case that the Hi can be regarded as self-adjoint operators on Hilbert space (then called 
“partial Hamilonians”), it was shown in |24| that the following consistency condition (first 
mentioned in m) is necessary and sufficient for the existence of a unique solution on R^-^ 
for the initial value problem ([6]): 




0 yj,k. 


( 8 ) 


^Note that for N > 2, this is a slight (and relativistically natural) generalization of what the Lorentz 
transformations of simultaneous configurations yield. 
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On the more natural domain =5^, which has a non-empty boundary, boundary conditions 
may be required for the uniqueness of the solution. However, this depends sensitively both 
on the operators Hi and on the space-time dimension [2411271128] . In the case that the Hi 
cannot be regarded as operators on Hilbert space, e.g. because they contain multiplication 
or differential operators with respect to time, different techniques have to be applied to 
prove the existence and uniqueness of a solution [281129] . 


2.2 Currents, probability density and Hilbert spaces 

In order for the multi-time equations to make contact to statistical predictions and to the 
formalism of non-relativistic QM, it is essential that ([5|) implies the existence of a positive- 
definite tensor current {xi, ...,xn) as a bilinear function of the wave function. The 

current is supposed to be conserved as a consequence of the wave equations, i.e. 

= 0, k = l,...,N, (9) 

where is shorthand for 

ax^ 

In free multi-time A^-particle Dirac theory, i.e. for 

Hk = ■= -nhl^kj + ruk'Jk (10) 


one has: 

jl^l-k-Ni^Xi,...,XN) = 1p{xi,...,XNhi^ ■ ■ ( 11 ) 

where ip = ■ ■ ■ 7 ^. 

The positive component of then is i.e. the usual density. Experimental 

predictions for an equal-time hypersurface are based on this density. More generally, 
the adequate probability density ps for a space-like hypersurface with normal covector 
field n{x) is expected to be [20] : 

Pt.{xi,...,xn) = j'^'^"'^^{xi,.:,XN)n^^{xi)---n^j^{xN)- (12) 

We emphasize that the existence of the tensor current is the central requirement for 
a multi-time theory to make sense for physics. Even though it may seem only a minimal 
operational requirement, it is in fact sufficient to formulate a full realistic relativistic quan¬ 
tum theory [181120] . 

With this picture in mind, one can dehne natural Wparticle Hilbert spaces associated with 
a space-like hypersurface S [dOlldl) : 

:= (L2(S^)®C^ (•,-)e). (13) 

An appropriate scalar product (•, •)s can be defined if there exists a bilinear form of two 
wave functions that generalizes the tensor current j^, i.e. if there exists 


j^^-^^['ipi,ip 2 ]ixi, ...,xn) with dk,i,J^^-^'^-^^[ipi,^p 2 ]{xi, ...,xn) 
for all solutions of ([S]) , and with 


0, k = 1, ...,N 
(14) 


( 15 ) 

For free multi-time Dirac theory, we have: 


f^'"^^['ipi,'p2]ixi,...,XN) = V'i(2;i,"-,a;Ar)7r •••7iv^V’2(a:i,...,XAr). (16) 
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Then the natural choic^ for the scalar product is given by: 

One can easily verify that for j[V’i;V' 2 ] given by (11611 the scalar product reduces to the 
familiar expression {(j), x) = / d^xi ■ ■ ■ d^XN <^^X for an equal-time hypersurface 
Furthermore, the scalar product has physical meaning by its connection with ps from eq. 
(fT^ . This can be seen as follows: Define 

IIV'lls := (18) 

Let A C and let 1a denote the indicator function of the set A. If HV’Hs = Ij then by 
eq. ([12]), ||lyi'i/>|||. gives the probability for a spatio-temporal conhguration of N particles 
on S to be in A m- 

We note that a similar expression for the scalar product has been suggested by Rizov, 
Sazdjian and Torodorov |251I32| . They define: 

i4>,x) ■= / da{xi)--- / da{xN) ■■■,XN)n^^{xi) ■ ■ ■ nf,^{xN)- (19) 

The crucial difference to eq. (HZD is that Rizov et al. allow the space-like hypersurfaces 
to be different. This is problematic because according to the argument at the beginning 
of this section, a multi-time wave function is naturally only defined on the set of space¬ 
like conhgurations. However, an element of Si X • • • X is in general not a space-like 
configuration. 

Furthermore, there is no reason why then should be equal to unity. This is 

because a crossing probability ps can in general only be defined on sets of the form S'^. 
The physical reason for this is that a time-like conhguration may actually correspond to 
two points on the world-line of the same particle. Besides these physical arguments, the 
mathematical consistency of the multi-time equations dHj) on the whole of may lead 
to very restrictive conditions on the possible operators Hi, e.g. excluding a multi-time 
formulation of QED models which (disregarding UV divergencies) is possible on =5^ |151I161 
[29] . The domain also plays an important role for the possibility to rigorously introduce 
relativistic interactions by boundary conditions at the space-time points of coincidence in 
a simple 1 -|- 1-dimensional model [281133] . 

Returning to the Hilbert spaces , we can understand a multi-time system ([5]) on =5^ to 
dehne a unitary evolution between different space-like hypersurface S, S' as follows |29][3I]: 
Denote by the restriction of a solution of ([5]) to S. Then 

^ V’R ^ ( 20 ) 

dehnes a map with HV'l^lls = ||s' = ||s' [M]. This way of formulating the 

time evolution puts the multi-time system of equations in the foreground and yields the 
Hilbert space picture as a byproduct. Note that reversing the train of thought and first 
dehning Hilbert spaces and the map does not in general yield a multi-time 

wave function on (which may be necessary for a consistent interpretation |181l21j L This 
can be seen as follows: Pick two space-like hypersurfaces S 7 ^ S' with S n S' 7 ^ 0. Then 
the Hilbert spaces and are essentially different and for 4> S we may have 

^Note that the natural range of integration is (Ei x . •. x Ejv) D expressing that the particles are 
always within the domain. However, as discussed in |28| . for d > 1 space dimensions the difference to (O 
is only a zero-measure set. 
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(Piq) + (</>)(?)■ 

The connection to the Hilbert space setting for the Schrodinger equation is given by re¬ 
stricting to equal-time hypersurfaces Sj in a distinguished frame and identifying all Hilbert 
spaces for different t |29) . This is possible for flat hypersurfaces without changing the 
scalar product. Also, because of this identification, the above-mentioned problem cannot 
occur because trivially there do not exist t with n Sf/ 7 ^ 0 . 

2.3 A no-go theorem for potentials 

For classical mechanics, the “no interaction theorem” seems to rule out a sensible physical 
theory of N directly interacting particles in the Hamiltonian framework [221123) . Pre¬ 
supposing canonical quantization rules, this may also have significance for a relativistic 
quantum theory. There do, however, remain doubts if this connection is robust enough 
or if a quantum theory does not simply offer further possibilities (see |34| for a discus¬ 
sion). Considering the Dirac equation, one may especially think of spin which does not 
have a classical analog. In view of the severity of a no-go theorem, it is important that 
its assumptions exactly fit the appropriate situation. A no-go theorem appropriate for 
multi-time Dirac equations was recently proven in |24) : 

Let the operators in eq. ([5]) be of the form 

Hk = + ( 21 ) 

where is defined in eq. cni) and Wk is a matrix-valued function of its arguments 

which only acts on the spin index of the A:-th particl^. Then the consistency conditions 
()5|1 are satished if and only if the multi-time equations are gauge-equivalent to the case 
Wk = Wk{xk), i.e. of purely external potentials. 

It is therefore an interesting challenge to find different mechanisms for relativistic interac¬ 
tions. Multi-time QFT models provide further possibilities for interaction terms by using 
creation and annihilation operators [laiinDITlEg] but they encounter the problem of UV 
divergencies. For 1 -|- 1 dimensions, one can construct a rigorous, interacting model by 
boundary conditions on the set of coincidence points of two particles in space-time |281I33| . 
The 2BD equations offer a third possibility by allowing the total momentum operator in 
the “potential” terms. This departs from the form ([5|) of the multi-time equations in that 
time derivatives are present in the interaction terms. Consequently, the system is not 
“Hamiltonian” and thus the “no interaction theorem” cannot be applied. The option of 
the 2BD equations seems especially attractive because the construction works for 1 -|- 3 
dimensions and does not lead to any known UV divergencies. This constitutes a strong 
motivation to analyze if these equations can indeed be used as the basis of a self-contained 
(and in this sense fundamental) relativistic multi-time theory for two particles. 

3 The Two-Body Dirac equations as multi-time evolution 
equations 

In this section, we introduce the idea underlying the 2BD equations. Starting from a gen¬ 
eral form of the multi-time equations which admits differential operators of any order in 
the interaction terms, we analyze the implications of a necessary compatibility condition 
when combined with Poincare invariance (sec. [321). The central features are that the in¬ 
teraction terms are related in a specific way and that a certain variable x_\_ appears which 

^This assumption is needed for the free Dirac current (HD) to be conserved in the presence of the 
potentials. 
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generalizes the spatial distance in the center of momentum frame in a subtle and covariant 
way. 

For definiteness, we then specialize to a general class of the 2BD equations for fermion- 
antifermion systems which was first presented by Sazdjian in [1] (sec. 13.31) . The focus is 
on the case of Dirac (spin-^) particles because for Klein-Gordon particles already in the 
one-particle case the candidate probability density can be negative |35| . The choice of 
the particular class of 2BD equations is made because (a) its compatibility can be seen 
straightforwardly, (b) it can in some sense be derived from the Bethe-Salpeter equation |36) 
and (c) there exists an analysis of its tensor currents [25]. We emphasize that the focus on 
the this class of the 2BD equations is not very restrictive as there exists a formal relation 
of Sazdjian’s 2BD equations and those proposed by Crater and Van Alstine [HIET]. The 
main features of these two classes of 2BD equations are shared, such as a dependence of 
the potential terms on the total momentum operator. 

Having introduced the equations, we comment on the question of how one can understand 
them as evolution equations for a multi-time wave function fsec. 13.5^ . As these equations 
leave the standard framework for the existence and uniqueness theory of quantum mechan¬ 
ical wave equations in several respects, it is required to introduce at least a preliminary 
mathematical understanding. 

3.1 Notation 

We use the following abbreviations: 

X = Xi-X 2 , X = {xi+X2)/2, 

d 

= A: = 1,2, 

P= {Pi-P2)/2, P=pi+p2. (22) 

3.2 Implications of the consistency condition for the general form of the 
equations 

Consider the following general form of multi-time wave equations: 

Diil;{xi,X2) = 0, i = l,2. (23) 

The crucial point is that in contrast to the form in eq. ([5|) (which led to the no-go theorem 
m) we now allow the operators Di to be differential operators of any order, including 
infinity. 

Recall from sec. 12.11 that certain compatibility conditions have to be satisfied in order for 
the wave equations (|23l) to have solutions. However, because the Di are not necessarily 
first-order differential operators, eqs. (|23p cannot be cast into the form Q. Therefore, the 
consistency condition ([5|) is not appropriate. For general operators Di, one cannot expect 
to hnd a simple replacement of eq. ([ 8 ]) which is also both necessary and sufficient. However, 
a necessary condition reads 

[Di,D2] = AiDi-I-A 2 D 2 (24) 

where at least formally the Aj can be operators. 

To see that (I24p is a necessary condition, consider a solution 'll: of (l23l) . Assume that one 
cannot write [Di,D 2 ] in the form ((241) . In general, we then have [Di, D 2 ] = X 1 D 1 +X 2 D 2 +R 
where R'ljj 7 ^ 0. Thus: [Di, D 2 ]'ip = RpJ 7 ^ 0 in contradiction with the fact that [Di,D 2 ]'ip = 
0 trivially holds for any solution of (I35p on which the action of D 1 D 2 and D 2 D 1 is well- 
defined. 



Remark: Note that in the case that the operators Di are of first order, condition (|24p 
bears some similarity with the consistency condition of Frobenius’ theorem from differential 
topology. However, it was discussed in |241 sec. 2.4] that the conditions for multi-time wave 
functions are, even for first-order operators, different from the statement of Frobenius’ 
theorem, one of the reasons being the number of components of ' 0 . 

Furthermore, note that ([241) applied to the first-order multi-time equations ([5|) seems to 
lead to a weaker condition than the previous condition ([8]). For (I24p it is sufficient for the 
right hand side to be a linear combination of the operators - Hk) instead of having 

to vanish. At first glance, the result in |24| . according to which ([ 8 ]) is a necessary and 
sufhcient condition thus seems to deem f[24p a too weak condition. However, |24| is tied to 
the case that the operators Hj. are operators on Hilbert space (or more generally operator¬ 
valued functions of the time variables). As the operators Di contain time derivatives and 
therefore cannot be regarded as operators on Hilbert space, it does not follow from |24j 
that the right hand side of eq. ([241) has to vanish. Still, the question of a necessary and 
sufficient condition for the system of equations ([23l) remains open. 

The important question now is: Which operators Di satisfy condition (1241) ? In order to 
reduce this question to a tractable form, we assume with |4] that operators Di take the 
following form: 

Di = A,o + f = l, 2 , (25) 

where Di^ are the operators of the free equation and V is an operator the structure 
of which is yet to be determined. Generally, it may depend on x, p and P as well as the 
gamma matrices qf , 72 • A dependence on X is excluded if one aims at a Poincare invariant 
theory. Note that (|25p introduces a relation between the interaction terms in the two wave 
equations ([2cip . We remark that the form (1251) may imply that the wave function of eq. 
([23p cannot be directly identified with the wave function of, say, the Breit equation, or the 
class of 2BD equations of Crater and Van Alstine. A wave function transformation may 
be necessary to relate the two types of wave functions (see mm)- Notwithstanding, we 
analyze the theory resulting from (I25p on its own terms. 

Using the form (I25p . we obtain: 

[Di,D 2] = DifiD2fi + dIqV + D2fiV D2fi + D2fiV DifiV 

— A> 2 ,o— D2,o^ ~ DifiVDi^o — DiflVD2fiV. (26) 

As the Di^Q are supposed to correspond to the operators in a free wave equation (acting 
only on the coordinates and spin indices of the z-th particle), we have [Di^, T* 2 ,o] = 0 and 
the first summands in the lines of eq. (1261) cancel. 

Aiming to bring (I26p into the form of the right hand side of eq. (I24D . we calculate the 
expression 


-[Di^o,V]Di + [D2fl,V]D2 = - DiflVDi^o - DiflVD2,oV + VDIq 

+ D2fiV £> 2,0 + D2fiV Di^qV — V D 2 Q. (27) 

Comparing eqs. (1261) and (EH, condition (I24p is satisfied if 

[Dio-Dio, V] = 0. (28) 

Specializing to the Dirac case, we evaluate (1281) for 

Difi = -fi-pi-mi. (29) 
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(30) 


Then: D^q = pf + mf and 


eq. ( 1281 ) reduces to: 

\Pi-P2:V] = 0 . 

Rewriting this equation using total and relative momentum operators yields: 

[P-p,V] = 0 . 


(31) 


Now, because of Poincare invariance, V must not depend on X. Thus, we arrive at the 
condition 




(32) 


To further evaluate eq. (|32l) . note that the only Poincare-invariant quantities involving x 
are of the form x ■ q where g is a quantity transforming as a 4-vector (e.g. x, P,p, 71 , 72 ). 
Thus, oc Qfj^ and (|32p requires 5 T P in the Minkowski sense. This can in general only 
be achieved if has the form 

(33) 

where ■= (1 - p~p') operator version of) the projection operator on the subspace 

“orthogonal to P” and q again transforms as a vector. 

The most important consequence of this can be seen by considering the case that q = x. 
Then eq. (I33p implies that the only dependence of R on x may be via 


x^, := 


(34) 


To see this in detail, note that x ■ q = x ■ x_\_ = x ■ (vrx) = x ■ (vr^x) = xx • 


On the meaning of xx: In classical mechanics, the analog xx of the operator xx 
(i.e. where P in eq. (|22p is replaced by a time-like 4-vector P) has the following meaning. 
Consider the relative spatial coordinate x = xi — X 2 between the two particles in the center 
of momentum (c.m.) frame, i.e. the frame where the total momentum 4-vector takes the 
form (P*^, 0,0,0). Then xx = (0, x). In this way, one can see that xx is the covariantization 
of ( 0 , x). 

Note that generalizing a non-relativistic law by the replacement (0,x) ^ xx would seem 
suspicious. If P were the total momentum of the total system, it would be inacceptable 
because in this context total quantities do not have any physical meaning but are only 
used to define coordinate systems. However, the use of xx is restricted to autonomous two- 
particle system^ which are to be thought of as subsystems of a larger system. Then, the 
total momentum P is meaningful. Yet, one might object against the use of any preferred 
frame, even if it is dynamically preferred, such as the c.m. frame here. However, this 
criticism is alleviated by the fact that the replacement ( 0 , x) —>• xx is never used in the 
derivation of the necessity of Xx- Rather, the crucial starting point is the for m (|25p of the 
operators Pj - which is far from directly assuming a covariantization of a non-relativistic 
law of motiorll. 


'^See |31| for an analysis of when an autonomous subsystem description is possible for relativistic systems 
with spin. 

®A related subtlety of the notion of Lorentz invariance was critically discussed by Bell in [38]. 
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Further remarks: 


1 . One can also motivate the necessity of x± very concisely in the context of two-body 
Klein-Gordon equations of the form (j5? — rnf — V)'il> = 0, i = 1,2 where F is a scalar 
and Poincare invariant potential (operator) O sec. II]. However, this derivation is 
of limited significance for the approach to the 2BD equations taken here because the 
square of the operators (I25p does not in general yield pf — mf — V with a scalar 
potential I^. 

2. The necessity of x± does not follow from the connection of the 2BD equations with 
the Bethe-Salpeter equation (see |36jl. Rather, the insight that x± is necessary to 
formulate differential equations of the type (I23p is itself used to make the so-called 
“relativistic instantaneous approximation” which creates a manifest xj_-dependence 
of the potential terms. 


3.3 Two-Body Dirac equations for fermion-antifermion systems 


For the rest of the paper, we now specialize to an important class of 2BD equations which 
is almost identical to one discussed above: the case of spin-^ particle-antiparticle pairs. 
This case is particularly relevant as one aims at a theoretical description of mesons and 
their spectra IMS]. Equations for fermions with the same charge are, on the other hand, 
not believed to describe bound states (as are the equations below) and therefore not to 
lead to particularly interesting subsystem dynamics. 

The class of 2BD equations for particle-antiparticle pairs (first introduced by Sazdjian 
in [H sec. VI]) is given by 


Di'il}{xi,X 2 ) = 7 i-pi - mi - (- 72 -^2 + m-2)F ' 4 ){xi,X 2 ) = 0 , 
D2ip{xi,X2) = 72-^2+ m2-b (71-Pi+ mi)F '0(xi,X2) = 0 . ( 35 ) 


Here, ip is a 16-component wave function. According to standard sign conventions, particle 
2 is the anti-particle. 

The form of the equations is motivated similarly to the approach via eqs. (|23l) and (1251) . the 
only difference being that one has to account for the symmetries of the fermion-antifermion 
system. More precisely, the first of the equations (]35|) has to be obtained from the second 
via charge conjugation and mass exchange [Up. 3411]. This changes some signs when 
compared with the form of the operators Di in eq. (1251) . 

V is an operator which may depend on P,p, 3 :_l and the 7 -matrices in a Poincare invariant 
manner. The symmetry of the fermion-antifermion system demands jH p. 3411]: 

V"(1,2;7i, 72) = F(2, l;-72,-7i) (36) 


where “1 0 2” indicates the exchange of particle labels in quantities like pj, Xi (but not the 
7 -matrices). We remark that here and in the following the notation V{...) is only meant 
to emphasize possible dependencies on certain variables. V is always the same operator. 
Then the 2BD equations (I35D satisfy also the compatibility condition (I24p because the 
following relation holds p. 3412]: 

[111,112] = - ill ■ Pi,V]Di + [^2 ■ P2,V]D2. (37) 

Special choices of V may yield [Zli,Il 2 ] = 0 jU sec. VH]. 

®Note that starting from a different point, Crater and Van Alstine were in fact able to derive 2BD 
equations with scalar interactions as “square roots” of corresponding scalar interacting two-body Klein- 
Gordon equations [3]. 
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3.4 Basic mathematical questions 

Due to the dependence of V on P, the 2BD equations (j35l) are of infinite order both in space 
and time coordinates. This immediately raises difficult mathematical questions which are 
generally not addressed in the literature (compare e.g. [3l[3]). Besides the question of the 
compatibility of the 2BD equations which was already discussed in sec. 13.21 one may ask: 

1. What are appropriate initial data? 

2. What is an adequate space of solutions? 

Concerning 1., note that in contrast to a first-order multi-time system ([5]), one does not 
expect initial data to consist only of prescribing the wave function for configurations on 
a space-like hypersurface like in ([6]). For wave equations of n-th order one would rather 
expect that also (n— l)-th time derivatives have to be prescribed. If this analogy extended 
to infinite order, this understanding of time evolution for the 2BD equations could not 
make sense since prescribing all derivatives of an (analytic) function on a Cauchy surface 
is equivalent to writing down the solution on the whole of K® (for two particles). One may, 
however, hope that the fact that the infinite order only arises from the dependence of V on 
the total momentum operator may help to identify sensible initial data (see section [33|) . 
With respect to point 2, note that the possible spaces of initial data are usually a good 
starting point for defining Hilbert spaces on which at least the non-relativistic existence and 
uniqueness theory is usually based. Due to the occurrence of powers of P to infinite order, 
and therefore of time derivatives, it is clear that this setting cannot be used without mayor 
changes. Moreover, the considerations in section [22] show that the natural Hilbert spaces 
depend via the scalar product (1171) (and corresponding statements about self-adjointness 
etc.) on the form of the conserved tensor current of the theory ([Uj). The question of 
conserved currents for the 2BD theory will be addressed in section HI 


3.5 A preliminary mathematical understanding 


In this subsection, we propose a way how one can understand the 2BD equations in a 
preliminary way for superpositions of eigenfunctions of the total momentum operator. 
Assume that the only momentum dependence of V is via P (explicitly via P^ or implicitly 
via x±). We write: V = V{x±, P). Let ifip be an eigenfunction of P, i.e. a function of the 
form 

'ippixi,X2) = 'ip{x)e~"^'^. (38) 


Then 

V{xp,P)il:p = H(xx,P)V’P, (39) 


where V{x_\_, P) is the matrix-valued function which is obtained by replacing P in x± with 
its eigenvalue P. In this way, we can regard H as a multiplication operator. 

For ipp, eqs. (1351) then constitute a first order system of differential equations, analogous 
to ([5D. The analogy with ([^ then suggests that adequate initial data are of the form of 
prescribing ijjp{xi,X 2 ) on a space-like hypersurface, e.g. for Xi=X 2 = 0: 


V’p(0 ,xi,0,X2 ) = '0o(x)e*^'^, xeM®. (40) 

The role of the two eqs. (I35p then is to (a) time-evolve ip in x^ and (b) determine P^. 
More generally, one should consider superpositions of eigenfunctions of P. These functions 
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are necessary to describe localized wave packet^ on the configuration space of two particles. 
Let Z denote further quantities, e.g. the relative momentum eigenvalues, which classify a 
suitable space of “relative coordinate wave functions” V'(x). Then: 

/ d"Pc(P,Z)^z(x)e-'^-^, (41) 

where it is understood that each is determined by eqs. (f35l) by demanding that 
be a solution for every Z, P. 

The further strategy of the paper is the following: Setting aside all further mathematical 
question^, we assume that solutions of the form (I4ip exist, at least for superpositions of 
finitely many eigenvalues of P. This permits us to analyze the central physical question if 
there exist adequate conserved currents for the 2BD equations. 

4 The question of current conservation 

Recall the central place of the tensor current [ijji,'4’2] and especially in the 

general structure of a multi-time theory (sec. 12.2p . The importance of j has also been 
recognized by various authors in the context of the Two-Body Dirac equations, in particular 
for the question how to construct scalar products and corresponding Hilbert spaces, see 
|321I39) for the spin-less Klein-Gordon case and [4l[25] for the Dirac case with spin. 

Here, we first review previous results for the 2BD equations, adding details and clearly 
stating critical assumptions (sec. KM- It turns out that the free Dirac current is not 
conserved and that one has to construct suitable replacements. While replacements can be 
found, they are neither unique nor simple. We follow Sazdjian |25) to pick a particular one 
in order to be able to further analyze the resulting theory. The expression for the scalar 
product (see eq. (|17p ) is in general not positive definite and a comprehensive analysis of the 
associated problems has (to our best knowledge) not been performed yet. Therefore, we 
discuss whether the “scalar product” can be rendered positive definite by (a) restrictions 
on the function space Isec. 14.2P and (b) restrictions on the potential terms Isec. 14.3P . The 
further implications of the respective changes are illustra ted in both cases. 

4.1 Previous results 

In the following we assume that V satisfies the following hermiticity condition: 

VHxl,P) = lh>ix±,Phh2 (42) 

where V (without the hat) was introduced in eq. (f3^ . 

Claim 1: Let ipp, ipp/ he eigenfunctions of P. Then the free Dirac tensor current 'ipp'] 

ippJi'y^f'ipp' is conserved if and only if no interaction terms —(—72 • P 2 + 'm 2 )V and 
(71 'Pi + mi)V, respectively, are present in the 2BD equations. 

^Note that the issue of localization is not as problematic for the 2BD equations as e.g. in relativistic 
quantum field theory where the Hamiltonian is assumed to be bounded from below. The reason is that, 
as with the single-particle Dirac equation, for the 2BD equations negative eigenvalues of the energy- 
momentum operators are possible. 

®A good starting for the question of existence and uniqueness might be to first specialize on the case of 
a total momentum eigenfunction (I38II . One may then hope that via m a suitable space of solutions can 
be constructed. 
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Proof: Consider 


(V>p7f72V’P') = - (7fPi,MV’p)'^7i7272V’P' + '0 p72(7i‘pi,mV’p')- (43) 

Denote^ Vipp by Vpif^p - which still contains a x_L-dependence. The first of the 2BD 
equations ([35]) yields: 


7i-pi'0p/ = [nil + {-^2 ■ P2 + m2)Vpi]'ippi. 

Using the relations (|42p as well as ( 7 (()^ = 7 fc 7 fc 7 fc h follows that 
{'n-piipp)^ = V’p 


Vp{m2 + 72 - m) + rni 


7i72 


(44) 


(45) 


where the arrow indicates the direction in which the derivative acts. Combining eqs. (1431) . 
fl44|) and (|45)), we obtain: 


idi,i,{'ipp'yi'y2'^p') 


'^p 


i — 

-Vp{m 2 + 72- p 2 )l 2 + 72 (-72 • P 2 + m 2 )Vp' 


P’P'- 


(46) 


We note the following points: 

1. The term with mi has dropped out. 

2. The term with m 2 , i.e. —Vpm 2 'y 2 + 72 ^^ 24 ^' yields zero only in the case that V does 
not contain 72 -matrices and for P = P'. 

i — 

3. Even in the latter case, the remaining term —Vp^ 2 ' P 2 72 — 72 72 ' P 2 Vp' does not 
vanish because [ 73,72 • P 2 ] 7 ^ 0. 

4. If P 7 ^ P' and V is not constant, not even special choices for V make the rhs. of 
(|46p vanish. The appearance of both Vp and Vp/ is unavoidable because the basic 
mechanism which allows the 2BD equations to circumvent the no-go theorems is the 
use of these momentum-dependent terms. 

5. The only case in which the free Dirac current is conserved is that the 2BD equations 
do not contain the interaction terms from the very beginning. 

Analogous reasoning for —id 2 ,u {P’plil^P’P') establishes the claim. | 

Claim 2 (see |25l p.l625]): There exist currents ji(,t[V’P; V’P'j which are conserved by the 

2BD equations. 


Proof: One looks for a term j^’^dl'^P^'^P'] such that 


Ant [ V'P, V’P' ] : = jZe [V’P, V'P' ] + Jadd I V'P, V’P' ] 


(47) 


is conserved. We leave away the square brackets [V'P,V’P'] in the following for notational 
ease but emphasize that in order to be able to treat U as a matrix this understanding is 
crucial. Defining 


Fi ■■= dipj^ 


piu 


F^-.= 


d 

dxi. 


•.'free ’ 


F := 81^^82,uji 


{lU 


®This replaces the previous notation V{x±, P) to fit in the equations. 


(48) 
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we see that has to be a solution of the equations 

= -Fr, = --ft- (■<£•) 

Such a solution is easy to find |25l eq. (3.19)]. Let G{x — x') be a Green’s functions of the 
four-dimensional wave equation, i.e. 

\JixG{x — x') = 6 ^^^{x — x'). (50) 

Then for any pair of Green’s functions Gi, i = 1, 2, a solution of (I49p is given by: 

Jadd(^i> ^ 2 ) : = -^1 J “ ^ 1)^1 ^ 2 ) -^2 j G2 {x2 - X2)F^ (xi, X 2 ) 

+ 8^82 j d^x'id‘^x' 2 Gi{xi — x'i)G 2 {x 2 — x' 2 )F{xi^x' 2 ). I (51) 

Remark: 1. From eq. (fMI) it is obvious that is not defined uniquely. One has to 

make a choice of the Green’s functions Gi, i = 1,2. Sazdjian’s choice i^^ Gi = 

Ga, i = 1,2 where Ga is the advanced Green’s function, with the reason that this 
would be “the only solution of (14911 which vanishes when the interaction is switched 
off” [251 P- 1625]. 

2. The construction of jint is very general (one might even say too general) since it 
would have worked for any F^^F^ defined by eq. ()48l) . So what is the significance of 
the existence of conserved jmt’s? A good answer would be to point out, for example, 
a unique current with the required properties like a positive component (see sec. l2.2p . 
However, a general argument why any of the possible dehnitions jint should yield a 
positive dehnite current simply does not exist. 

3. Nevertheless, one may ask the question whether given further assumptions the cur¬ 
rents are positive dehnite. Further assumptions might even be plausible, for example 
if they concern special potentials V. In the end, it is only important that the cur¬ 
rents are positive dehnite for realistic choices of V. An approach involving further 
assumptions was chosen by Sazdjian |25| which we shall critically review next. 


Sazdjian’s paper does not directly deal with the question whether the currents jint are 
positive dehnite but rather with the one whether the associated scalar product (see eqs. 
czD, (HU)) is positive dehnite. However, these two questions are equivalent as long as 
the wave functions admitted in the construction of the tensor currents (and in the scalar 
product) are not subject to restrictions which forbid localized wave packets. 

Sazdjian’s results for the scalar product derived from according to (|17p are as follows. 
The above-mentioned choice of Green’s functions leads to the following expressioiJ^ for a 
scalar product, for two eigenfunction^^ of P and the special case of S = i.e. 

an equal-time hypersurface with normal covector held n = (1,0,0,0) |25[ eq. (5.11)]: 


{pP,Pp')T.t 


L 


da{xi)da{x2) ji^P'ipp,'iljp']ixi, X2)nf,{xi)nu{x2) 


StxEt 

lini f fip - Hp,7i72Vp 

e^O I 


+ (P°' + Pp 


,0n yp'+ien - bp —i 


po' -pop 2 i£ 


fit 


(52) 


^^Sazdjian only uses eq. USD) with Gi = G2. 

^^We have adopted our notation conventions. Besides, the range of integration is corrected according to 
the remark in sec. [U below eq. Gi) such that the integration over xi and X 2 is over the same equal-time 
hypersurface instead of two different ones. 

^■^For more general wave functions such as in (EH, the definition of (■, ■}Et can be extended by linearity. 
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where (•)* denotes complex conjugation (without transposition). The limit e —)• 0+ comes 
from the definition of Ga{x) by its Fourier transform 


Ga{x) 


lim 

e->-0+ 


d‘^k 

(27r)^ — 2ik^e 


(53) 


Let ipp = e (pi{x), V’P' = e (p 2 {x). Then for Sazdjian obtains from 

eq. (1521) EH eq. (5.12)]: 


(V’P,'0P')St 


(27r)35(3)(p_p') J d^^(j)^{x) 


7?72°-Vp7?72Vp + 4(p0)2 


dVp ' 
5(P2) 


Mx) 

(54) 


where x = (x^ = 0,x). We note that because of the delta function one can use P = P' 
everywhere inside the integral. 


Remark: Ignoring the delta function, for P = P', eq. (1541) should yield the square of a 
norm. It therefore has to be positive. However, in general (i.e. independently of Vp) only 
the first term in the square brackets of eq. O, which corresponds to the usual expression 
rW;t^ = f in the Dirac case, yields a positive contribution. Thus, as recognized 

by Sazdjian Ea p. 16310: 

If the potential V is explicitly independent of P^ in the c.m. frame, the expres¬ 
sion of the norm (1541) shows that its kernel still depends on V. This implies 
that V must satisfy some inequality conditions to guarantee the positivity of 
the norm. This question was examined in more detail in Ref. sec. VII A]. 

In Ref. El sec. VII A], it is suggested to make a wave function transformation El (7-1)] 
which would map the norm given by eq. (1541) to the free / norm. However, this is only 
possible for operators V’g for which Vp does not depend on P^ in the c.m. frame El P- 
3423] and if in addition the following condition El e<l- (7-6)] is satished: 


iTr(7i- 


P 


72 • 


P 


V) < 1. 


(55) 


Presumably, the trace is to be taken over the spin components of V. The question whether 
the independence of Vp of P^ in the c.m. frame is a reasonable assumption is not clarihed 
in El- Tbe more recent article [HI eqs. (A4), (A9)] even seems to show the contrary. 
However, set aside this confusing point, a much more basic question is left open. Because 
V = Vp and because P is a property of ijjp, i.e. of the wave function, it is unclear how one 
should regard conditions that lead to the positivity of (1541) : 

1. as conditions on the space of admissible wave functions, or 

2. as conditions on the operators V, given their domain? 

These questions are discussed in none of the references EllBlES] for the 2BD case. 

In the following subsections we analyze the consequences of these two possibilities (see sec. 
1221 for possibility 1 and sec. 14.31 for possibility 2). 

^^For clarity, notation and references in the quote are adapted to our conventions, without changes in 
content. 
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4.2 Can the positivity of norm and currents be guaranteed by restriction 
of the function space? 

Recall from secti on 12.21 that the hope is to be able to regard the 2BD equations as defining 
an evolution mapo for any pair of space-like hypersurfaces S, S', i.e. 

^ V-s ^ V's' (56) 

with = 'i’T.'iq) if g G S n S'. should be unitary in the scalar product defined 

by jint according to eq. m- 

( 2 ) 

However, as we saw above, this construction does not yield a scalar product on ' because 
it is in general not positive definite. Thus, we define 

"^1;“ := {(/> G :(()), ()))s < oo A ((/>,()))s > 0} U {0} (57) 

(< 2 \ 

as the subspaces of for which (•, •)s is actually positive-definite. The question is: does 
define an aeceptable space of functions? 

To decide on this question, consider the following points: 

1. It is not clear anymore that contains all physically relevant functions. One can 
see this e.g. from (|4ip and (1541) . Any reasonable quantum mechanical theory which 
describes matter should be able to describe localized wave packets. To construct 
these wave packets, one in general requires all Fourier modes tpp as in (|¥T]) . However, 
for a general V, the requirement of positivity of (I54D implies conditions on the P’s 
such that some are not admitted in the construction of wave packets. Furthermore, 
these conditions are mathematically quite involved and do not serve a clear physical 
purpose. 

2. The do not, in general, define Hilbert spaces. Completeness may be violated. 

Even worse, the may not even be vector spaces. Consequently, the mathemati¬ 
cal structures on which both the usual quantum formalism is built up, break down, 
including the self-adjoint operator observables as well as the standard approach to 
define the time evolution. Of course, one may consider to further reduce the admis¬ 
sible functions by replacing with some Hilbert space Td* contained in it. In 

fact, a similar route was suggested by Sazdjian [251 p. 1624]. This, however, further 
strengthens the criticism of point 1 and still leaves open the question if the usual 
mathematical structures can be defined. 

A simple analogy: To appreciate the problems that accompany the restricted function 
spaces Td^^ and Td* , consider the following example. Let the Hilbert space of our theory 
be given by with “scalar product” 

(u, w)a '■= Aw, where A = ^ J 

Of course, (•, ■)a does not define a scalar product on C^. So in analogy with (I57p we define 

:= {u G : uUu > 0} U {0} (59) 

'^^Note that this question is independent of the fact that the operators Di appearing in the 2BD equations 
are not operators on Hilbert space. 
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as the subset on which (•, ■)a actually is a scalar product. We note that e.g. (0,1) (which 
may be a physically relevant vector to represent a spin state) is not contained in (cf. 
point 1 ). 

Moreover, is not a vector space, because the for v\ = (1, G and V 2 = 

(—1, the sum vi + V 2 = (0,1) is not an element of Furthermore, it is 

also not complete, as the following example illustrates: Consider the sequence given by 
Vn = (1,1 — 1/n). We have: {vn,Vn)A = 1 — (1 — 1/n)^ > 0 and thus Vn G However, 

the limit v = (1,1) has norm zero, i.e. v 0 These problems are analogous to point 2 

from above. Note that they can be overcome by dehning even smaller Hilbert spaces Ti* as 
the span of ( 1 , 0 ) instead. ("H* is a complete vector space for which (•, ■)a dehnes a scalar 
product). However, even more physically interesting vectors get lost this way. 

To extend the analogy, suppose that the “wave equation” of our theory is 


i—u = Bu 
dt 


where B is the following 2x2 matrix 


(60) 


B 


0 i 
-i 0 


(61) 


The matrix B is self-adjoint with respect to the canonical scalar product on but not 
with respect to (•, ■)a- Thus, it defines a time evolution on but not necessarily on 
Let us analyze the consequences. Given u{t = 0) = uq, we have 


u{t) = exp{—iBt)uo = (I 2 cost — iB smt)uo. (62) 

Let tio = (1,0) G n T-L*. Then u{t) = (cos t, — sin t) which is in general not an element 

of (neither of T-i*). 

One may try admitting only initial data uq = (a, b) G for which also u{t) G Vt. 
Then |a| > | 6 |. We have u{t) = (acost + 5sint, —asint + boost) and therefore 

u^{t)Au{t) = (|op — | 6 p)(cos^ t — sin^ t)-I-4Re ( 0 * 6 ) cost sint. (63) 

We ask: Do a,b G C with |a| > | 6 | exist which make this expression positive for every 
t? For an answer, consider (1631) for (i) t = 7 r /4 and (ii) t = Sir/A. In case (i), we 
have sint = cost = lly/2 and obtain as a condition that Re (a* 6 ) > 0. In case (ii), 
sint = — cost = l/\/2 and we obtain the condition Re(a* 6 ) < 0, in contradiction with 
(i). We conclude that the restriction to is not in any way consistent with the given 

time evolution ((601) (neither is the restriction to %*). This illustrates the problems with 
dehning the time evolution from point 2 above. 


Comparison with Klein-Gordon theory: If the above analogy extends to the case of 
the 2BD equations, the logical consequence is to reject restrictions on the function space. 
However, in view of previous claims about the consistency of a Hilbert space picture for 
interacting two-body Klein-Gordon (KG) equations in [2511321139) . one may wonder if points 
1,2 are actually as severe as they seem to be. 

These sources (especially [251 sec. HI]) show the following. In the KG case, one can identify 
potentials such that the scalar product given by the conserved tensor currents of the theory 
according to (|17p is positive dehnite on a subspace T-L* of the Hilbert space By. of the theory. 
Then B* is again a Hilbert space, corresponding to one of four possible choices of the sign 
of eigenvalues of the operators pi ■ P/P‘^ and p 2 ■ P/P^. One may thus hope that the 
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problems of point 2 do not appear. 

However, this approach seems to disregard problem 1: the subspace T-L* does not contain all 
physically relevant functions. One cannot, for example, represent localized wave packets by 
wave functions in %*. To do so would require basis vectors from all of TLy.- Furthermore, 
problems with the self-adjointness of operator observables may occur (see also [321 p. 66]). 
A completely analogous situation is encountered in free one-particle KG theory |35l chap. 
3]. In this case, one draws the logical consequence that the KG equation theory cannot be 
considered a self-contained one-particle theory. By the same arguments, one also has to 
reject the approach via Ti* towards interacting two-body KG theory. 


Conclusion: One may wonder whether or not the situation for KG theory has any sig¬ 
nificance for the 2BD theory. As remarked after the discussion of the meaning of x±_ in 
sec. 13.21 the square of a the 2BD equations does in general not yield an interacting KG 
equation. Moreover, recalling the quote at the end of section 14.11 the implications of the 
two-body KG theory on the 2BD theory are limited. For the 2BD case, the “scalar prod¬ 
uct” is not positive definite even if the norm is independent of in the c.m. frame (see 
also [H p. 1627-28, 1631]). Furthermore, we note that the approach in the KG case in¬ 
volves both restrictions on the potential operator V as well as restrictions on the function 
space. The restrictions on the function space turned out to be inacceptable whereas there 
is no reason to reject restrictions on the potentials as long as they include the ones used 
in applications. In view of this situation, together with points 1-2 (as illustrated by the 
analogy), we conclude that restrictions on the function space are also inacceptable for the 
2BD theorj{^. The question if, on the other hand, there exist sensible restrictions on the 
potentials such that the scalar product is positive definite is the subject of the next section. 


4.3 Do special operators V exist for which scalar product and currents 
are always positive definite? 


Consider eq. (|54p “in the c.m. frame”, i.e. for P = (P'^,0,0,0). Then x_l acts as the 
multiplication operator with the spatial part x of the relative coordinate. Demanding that 
('^P,V’p)st be positive for all eigenfunctions of P, we obtain the following condition 
for Vp: 


(/.(x) 


7?72°-Vp7?72Vp+4(p0)2 


dVp ' 

d(p2) 


(/)(x) 


> 0 . 


(64) 


This condition should be satisfied for a reasonably general class of functions (p, e.g. for all 
4 > G L2(]^3^ (g) (jgi6^ ^]^g following points: 


1. Vp, which depends on P via the quantities x± = (0,x) and P^ = (P*^)^ in the c.m. 
frame, has to be boundec0 with respect to (a) P^ and (b) = —x^. (a) is easy to 

achieve, e.g. in the case that V does not depend on (P®)^ in the c.m. frame, (b) is a 
real restriction. We shall see the consequences below. 


2. Do solutions Vp of eq. (1641) exist? To answer this question, consider the class of 
scalar functions Vp = f{—x\). In the c.m. frame they take the form /(x^) which is 
independent of (P^)^. Thus, the term 4 (pO) 2 ^^Fp^ gq (|g^ vanishes. Making use 
of the fact that /(—x^) is real-valued as a consequence of eq. (]42p in the scalar case, 

^®One may even hope that the situation is better in the 2BD theory in the sense that the free Dirac 
current is positive definite, as opposed to the free KG current. 

^®“Bounded” in the context means that the absolute values of the eigenvalues of Vp = V{x±,P) are 
bounded. 
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( 65 ) 


condition ([ 61 ]) reduces to: 

())^(x) [1 - |/(x)|2] (/)(x) > 0 V(/.G L2 (m3),^c^6^ VxGM^ 
^ |/(x2)| < IVxGM^ 


Thus, we conclude that there do indeed exist special operators, e.g. V = /(—with 
|/(y)| < 1 Vy G M, for which the scalar product is positive definite on a general function 
space, e.g. for superpositions of eigenfunctions of P (|4ip with suitable drop-off conditionJ^ 
for \x\ \ —)• oo. Given any smooth and real-valued function g{y), such a function / can be 
constructed as 

f{y) ■= tanh 5 (y). (66) 

One may, however, ask: are these restrictions on V physically reasonabl^^? 

We try to answer this question by comparison with realistic potentials derived from quan¬ 
tum field theory in | 8 l appendix A]. One such possibility for scalar interactions in lowest 
order perturbation theory is | 8 l eqs. (A4), (2.17), (2.20)]: 




tanh 


1 glg2"^P 
2v^ 47r 



(67) 


where gi,g 2 G M, /i > 0. The question is: does Vi fulfil the positivity condition (1641) ? 

We first note that (i.e. Vi where P is replaced by an eigenvalue P) does indeed 

explicitly depend on even for P = (P'^, 0,0,0). This feature is shared with other 

possible potentials derived from QFT IseJ^ (S] appendix A]). Thus, we can neither use 
the simplified condition (1651) and nor the before-mentioned condition (f55|) of Sazdjian. 

Let us evaluate condition (I64D for Vi^p for eigenfunctions of P and in the c.m. frame. Then: 

P^ = {P^p and 'sj—x\_ = |x|. We have: 


dVi,p 

d(p 0)2 


l|^0|-3 glg2e 

4 dvr |x| 


cosh^ 


1 


1 .<7l.<72 e 

2\P^\ 471 |x| 


Let 

_ 1 9192 

^ 2\P^\ dvr |x| 

Evidently, y > 0. Then eq. (I64p becomes 


( 68 ) 


(69) 




1(1 - tanh2(-y)) + 27^7^--f--- 

cosh {—y) 


> 0 Vy > 0, V()) G C^®. 


(70) 


Note that 7 ® has eigenvalues ±1 for each j. Thus, eq. (f70l) yields the two conditions 

1 — tanh^(—?/) ±-^- >0 Vy > 0. (71) 

^ ’ cosh 2 (-y) - " ^ ^ 

However, for the function h{y) := 1 — tanh^(—y) — is negative for y > ^, 

corresponding to |x|e^l^l < Consequently, wave functions with internal part ^(x) 

with support concentrated around |x| = 0 have negative “norm” and “probability density”. 


^^Note that because of the form of the kernel of the scalar product (O the drop-off conditions may 
become modified as compared to the case V = 0. If this should turn out problematic, one could easily 
avoid the situation by demanding that Vp falls off with \x^\ —>■ 00 sufficiently fast. 

^®Note that when comparing such a bounded V with, say, a Coulomb potential (which is unbounded), 
one may have to take into account a wave function transformation (see [8] and the remark below eq. dSSJ). 

^^Note that also for Crater and Van Alstine’s form of the equations the potentials explicitly depend on 
VP2, called w in the references (see e.g. 0 appendix A]). 
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Comparison with the norm used by Crater and Van Alstine: Building on Sazd- 
jian’s work |25| . Crater and Van Alstine also considered the question of an adequate 
norm [Up. 9]. Their derivation of the norm is based on the following wave function 
transformation between the wave function ijj of eq. (f35ll and the wave function xf: appearing 
in the so-called “hyperbolic form” of of their equations [3 eqs. (52), (53)]: 

xjj = cosh(A)'0 (72) 


where 

A = tanh“^(V). (73) 

We note that this transformation is not a simple mathematical object because it evidently 
depends on the operator V. We continue with the assumption that it does indeed exist 
(which may yield further conditions on the potentials or on the function space) and analyze 
the consequences for the norm. 

Employing the transformation ([72]) for eq. (l5l|) . Crater and Van Alstine obtain (see [7] 
and m appendix B]; the result is adapted to our notation): 


('0p,'0P')st = (2vr)' 


>h(3)(p-p') I 




1 - 


2 „, 0„,0 


5(P2) 


h{x) (74) 


where the (j)i are dehned analogously to the (pi in eq. (I54p and Ap is the operator A 
with P replaced by its eigenvalue P. The symbol (•) indicates that the wave function 
transformation (I72p has been made. 

Considering eq. (1741) . we note that the expression for the norm (i.e. for P = P') reduces to 
the usual f |'0|2-expression for energy-independent potentials and is then positive without 
further restriction on the potentials. However, as evident from both jS] appendix A] and |9l 
appendix A], realistic choices of the potentials explicitly depend on the energy '/P^. Thus, 
equivalent restrictions on the potentials as given by condition (I64p also appear following 
Crater’s and Van Alstine’s approach. This, of course, has to be the case if the wave function 
transformation is to yield an equivalence between the 2BD equations of Sazdjian ([35D and 
the 2BD equations of Crater and Van Alstine. 

More precisely, one can see from eq. (I74p that the following condition has to be satisfied 
by A: 

<^^(x) 

for all ^ € L2(M^) (g) 

We now evaluate this condition for a realistic choice of A, corresponding to Vi from above 
(see eq. dSZ]), 13 eq. (A4)]). Then: 


1 - 4P^7i"72 


2 ., 0„,0 


a(p2) 


(^(x) > 0 


(75) 


Ai = tanh ^(Vi) 


1 glg2"^P 
2\/^ dvr 



(76) 


After a short and elementary calculation similar to the one leading from eq. dMI) to eq. 
m, condition (|75p reduces to 


X e‘ 


/i|x| 


< 


glg2 

dvr 


which is the same condition as before, with the same consequences. 


(77) 
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Conclusion: It is in principle possible to guarantee the positive definiteness of the scalar 
product by special choices for the potential operator V. This is particularly easy to achieve 
in the scalar case and for y’s which are independent of . Realistic choices for V such as 
Vi from eq. ([671) . however, are not independent of P^. This has the consequence that Vi 
does not satisfy the condition (I64p for positive definiteness of the scalar product and of the 
probability density. One may suspect that other realistic choices for V might suffer from 
the same problem. Therefore, they might not lead to a self-contained quantum mechanical 
two-particle theory which can possibly make statistical predictions in its own right. 

5 Gauge invariance 

In this section, we briefly comment on how the unusual mathematical structure of the 2BD 
equations influences the notion of gauge invariance. 

According to the view put forward in sec. 12.21 one should regard the tensor current ip], 

not the wave function ip, as the physical object. Transformations ip which leave j invariant 
are considered pure gauge. In the case of free multi-time Dirac equations, j^''[ip,ip] = 
ip'j^'y^ip and the gauge transformations are given by (see also [24] 1: 

iP{xi,X2) ^ (78) 

In the 2BD case, however, the tensor currents in the 2BD case have to be modified (see 
sec. mi). The possible replacements are momentum-dependent, i.e. their form depends 
on the wave function itself. Consequently, the form of a gauge transformation changes. 
In particular, the standard transformations (fTHl) cannot in general be considered gauge 
transformations anymore, because e.g. ip i—)• e~^^'^ip may change the eigenvalue P of P 
and P in turn is crucial for the form of jint- The class of gauge transformations is thus 
reduced to transformations 

6 {xi,X 2 ) = 0{x) (79) 

which do not involve the coordinate X on which P acts. Note, however, that the general 
gauge transformations (1781) may introduce terms into the multi-time equations (I23p which 
are not Poincare invariant. This is not possible using only the restricted class of gauge 
transformations (EH). 

6 Discussion 

In this paper, we critically reviewed various questions around the 2BD equations, rectifying 
critical issues where needed and pointing out unsolved problems. First, we presented in 
detail how the 2BD equations connect to the important work of Dirac, Bloch, Tomonaga 
and Schwinger on the multi-time formalism, as well as to recent developments. At this, we 
presented a general framework how to understand multi-time equations both mathemati¬ 
cally and physically. The latter aspect led to the requirement of conserved tensor currents 
with a positive component which can play the role of a probability density. The importance 
of this requirement cannot be stressed enough because it constitutes the link between the 
mathematical formalism and the statistical outcomes of experiments. Moreover, we hinted 
at the connection to realistic relativistic quantum theories. 

Starting from the general framework, we showed how the 2BD equations avoid a recent 
no-go theorem which, contrary to older well-known results such as the so-called “no in¬ 
teraction theorem”, is directly applicable to multi-time Dirac equations, not only up to 
quantization. The relevant mechanism for this is the inclusion of arbitrary powers of the 
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total momentum operator in the interaction terms. This allows to satisfy a necessary com¬ 
patibility condition by assuming a certain relation between the interaction terms. The 
compatibility condition then implies the appearance of (an operator version of) the vari¬ 
able x±. In classical mechanics, x± is the covariantization of the spatial distance in the 
center of momentum frame. The subtle arguments why the use of such a covariantization 
can be considered in agreement with relativistic physics were carefully discussed. 

We continued by focusing on a certain, still fairly general class of the 2BD equations for 
particle-antiparticle pairs which was first suggested by Sazdjian. Here, we pointed out the 
most important mathematical questions (which have not recieved sufficient attention in the 
literature) and proposed a preliminary mathematical understanding of the 2BD equations 
for superpositions of total momentum eigenfunctions. The question of the existence and 
uniqueness of solutions of the 2BD equations, however, remains open. 

The main part of the paper was devoted to the above-mentioned critical question of whether 
conserved and positive definite tensor currents exist. While the free Dirac currents are not 
conserved, there do exist possible replacements. These replacements are not unique and 
there is no general argument why for any of them the currents should be positive definite 
for arbitrary potentials and on a general function space. This situation motivated the 
question whether the currents can be rendered positive definite by restricting the function 
space or the admitted class of potentials. Our analysis started out from a previous one by 
Sazdjian which was, however, shown incomplete. In particular, Sazdjian did not discuss the 
question if further restrictions to render the currents positive definite are to be regarded 
as restrictions of the function space or of the potentials. 

These two aspects were the topic of the next two sections. First, we showed in detail that 
restrictions of the function space are not acceptable. The reason is that, roughly speaking, 
Fourier modes are excluded which are necessary for the representation of localized wave 
packets, for the self-adjointness of operators and for the usual way of defining the time 
evolution. We also commented on existing literature on similar questions, finding that a 
comprehensive analysis has never been done before for the Dirac case (the Klein-Gordon 
case anyway not leading to positive-definite currents). 

Second, we analyzed the implications of the requirement of positivity of the currents on the 
allowed form of the potentials, given a sufficiently general function space. The results were 
twofold: on the one hand, we found that it is indeed possible to identify a general class of 
potentials with the desired property. On the other hand, potentials which were suggested 
as physically accurate in the literature may in fact violate the requirements for positive 
definite currents. It should be emphasized that in any case the form of the probability 
density changes as a consequence of the fact that the form of the currents depends on 
the chosen potential. This is of relevance also for applications, e.g. concerning transition 
rates. (For spectra, on the other hand, the form of the probability density is unimportant.) 

As there are several possible motivations to study the 2BD equations, the implications 
of the results on restrictions for the admissible potentials can be regarded in different 
ways. On the one hand, for applications where realistic potentials are required, doubts are 
raised that phenomenological calculations of meson spectra based on the 2BD equations 
like in [HiiHiiiaiiniiii] do have a theoretical justihcation. To resolve the doubts would 
require to check the positivity condition (I64p (or, equivalently, (1751) 1 for the potential used. 
This, however, is not done in the literature, nor does an awareness of the problem seem 
to exist. In fact, we found that for a simple choice of the potential which was derived by 
Sazdjian from the Bethe-Salpeter equation, it is violated. The fact that this topic has not 
received attention before is somewhat surprising, considering that one of Crater’s and Van 
Alstine’s main motivations to introduce the first version of the 2BD equations was that 
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the Bethe-Salpeter equation possesses negative-norm states - and therefore does not have 
a clear physical interpretation [3]. 

Furthermore, assuming that physically realistic potentials could be found which also satisfy 
the positivity condition, the modihed probability density as compared to the |^p-density 
seems unusual. It would therefore be interesting to subject the modihed density to ex¬ 
perimental tests, for example by determining transition rates. In this respect, Crater and 
coworkers have investigated decay rates of quarkonium and positronium into two ohotonJ^. 
considering the effects of the modihed norm ([731) |40j . The theoretical results obtained com¬ 
pare well with other phenomenological approaches, but still lie outside of the error bars of 
the experimental data in all cases. These differences between theory and experiment are 
particularly interesting, considering that they appear at a place which is critical from a 
purely theoretical point of view: as we stressed in sec. 14.11 the tensor current is not unique, 
requiring the choice of two Green’s functions in eq. m- Note that the theoretical results 
for the mesonic spectra given in m, which are independent of the exact form of the tensor 
current, fit much better with the experimental data. One could take these findings as a 
motivation to study the question whether modifying the potentials or (as is particularly 
interesting) making a different choice of the tensor current could improve the theoretical 
results. However, before immediately drawing the consequence that such modifications are 
required, one should not forget that further (possibly critical) assumptions are involved in 
the process of calculating decay rates via the 2BD equations. This is obvious from the fact 
that the 2BD equations as a strict two-particle theory do not, by themselves, accomodate 
processes with variable particle numbers. A theoretical justihcation to nevertheless calcu¬ 
late decay rates using solutions of the 2BD equations therefore cannot be contained in the 
framework of the 2BD equations alone but has to come e.g. from quantum field theory. 
On the other hand, for foundational aspects in relativistic quantum theory, it seems re¬ 
markable that there do exist interaction terms for multi-time equations at all which satisfy 
the minimal requirements of Lorentz invariance and compatibility with a probabilistic 
meaning of the wave function. One class of these equations is given 


{ 71 -pi - mi - (-72 • P 2 + W 2 ) tanh [^(-x^)] }'i/;(xi,X 2 ) = 0, 

{72-P2 + 'm2 + {'yi-pi + mi)tanh[g{-x\)]}ijj{xi,X2) = 0, (80) 

where g{y) is an arbitrary smooth and real-valued function. The expression for the associ¬ 
ated positive tensor current is rather lengthy and can be calculated via eqs. (SHD, dSI]). The 
corresponding scalar product, evaluated on equal-time hypersurfaces of a special frame, is 
given by (1541) . As stressed above, the tensor current involved in this construction is not 
unique. Such a non-uniqueness of the currents is, however, not an uncommon situation 
in quantum physis. One can for example always add a term which is divergence-free - 
and sometimes this is even appropriate. Moreover, the additional freedom in choosing a 
Green’s function in eq. ([5ip might (in more general situations than (|8UI) 1 help to reconcile 
experimental and theoretical results for decay rates. 

Finally, one may wonder whether a similar approach as for the 2BD equations can be taken 
also for N > 2 particles. However, appropriate wave equations in a closed form have never 
been found. This may be due to the fact that there does not exist a generalization of the 
variable x± for N particles which allows to satisfy the necessary compatibility condition 
of the wave equations in a similar way as for two particles [371142] . 

am grateful to H. W. Crater for pointing this out to me. 

■^^Note that a similarly looking class of 2BD equations was suggested by Crater and Van Alstine [71 eqs. 
(52), (53)]. Eq. (I8UI) is a subclass of these equations for which the positivity of the scalar product and 
currents has been checked in sec. For the general class in |3, positivity may be violated. 
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